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(S| ■ ABSTRACT 

The formation scenario for massive stars is still under discussion. To further constrain current theories, it is vital to spatially resolve 
the structures from which material accretes onto massive young stellar objects (MYSOs). Due to the small angular extent of MYSOs, 
one needs to overcome the limitations of conventional thermal infrared imaging, regarding spatial resolution, in order to get obser- 

l/~) ' vational access to the inner structure of these objects. We employed mid - infrared interferometry, using the MIDI instrument on the 

y—i , ESO/VLTI, to investigate the Kleinmann - Wright Object, a massive young stellar object previously identified as a Herbig Be star 

precursor. Dispersed visibility curves in the N-band (8- 13 ^/m) have been obtained at 5 interferometric baselines. We show that the 
mid - infrared emission region is resolved. A qualitative analysis of the data indicates a non - rotationally symmetric structure, e.g. 
the projection of an inclined disk. We employed extensive radiative transfer simulations based on spectral energy distribution fitting. 

rf) , Since SED - only fitting usually yields degenerate results, we first employed a statistical analysis of the parameters provided by the 

radiative transfer models. In addition, we compared the ten best -fitting self -consistent models to the interferometric observations. 

i-^ ■ Our analysis of the Kleinmann - Wright Object suggests the existence of a circumstellar disk of 0. 1 M Q at an intermediate inclination 

ft' of 76°, while an additional dusty envelope is not necessary for fitting the data. Furthermore, we demonstrate that the combination 

^ of IR interferometry with radiative transfer simulations has the potential to resolve ambiguities arising from the analysis of spectral 

» I , energy distributions alone. 

-*— > . 
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•^ 1 . Introduction On the other hand, analysis of the data obtained by interferomet- 

^^ ric techniques, specifically from the dispersed visibility curves, 

One of the central unsolved issues in the study of massive is far from trivia i. This holds in particular wnen the number of 

O ; star formation is to which extent the accretion mechanism is base lines is small. In this respect, Monte Carlo radiative trans- 

'nT , comparable to the better - understood processes in low -mass fer ( M CRT) simulations represent an important tool to analyze 

• star formation, e.g., accretion of material from a disk struc- visibility curve s, by calculating synthetic curves for comparison 

O ture by the forming star at certain evolutionary stages. Although with the mo d e led to the observed data. 
<^> , there are several theoretica l scenarios for massive star formation 

^H . dZinnecker & Yorkdl2007h . the lack of observational data on the We are currently conducting a large interferometric study of 

'<. ■ inner structure of massive young stellar objects (MYSOs) pre- MYSOs with the Very Large Telescope Interferometer (VLTI) 

• £h vents a deeper understanding of massive star formation. There (ILinz et al.||2008|). Most of the targets were taken from a list of 

X are two main reasons for this problem: first, these objects are BN-type objects (|Henning et al.l[l984. First results on the ob- 

»H ' usually deeply embedded, so observations need to be done in J ect M8E-IR have been reported by Linz et al.| (|2009|). In the 

! the mid - infrared (MIR) or at even longer wavelengths. Second present paper, we concentrate on the Kleinmann - Wright Object 

the angular extent of these objects is rather small. In order to (|Kleinmann & Wright] |1973| - also known as M17 SW IRS 1 

resolve the inner structures of YSOs in general, high angular hereafter KWO). Based on NIR spectroscopic data, Porter et al.| 

resolution observations are required. Since diffraction - limited <H228|) classified the KWO as a transition object between a 

observations, even with 8 m-class telescopes, result in resolvable deeply embedded MYSO and a Herbig Be star. |Chini et al.| 

angular scales of a few hundred milliarcseconds, interferometric ( 2004 ) presented for the first time iJHK images, in which the 

techniques need to be applied to assess the innermost structures. KWO ls clearl y resolved into two separate sources. They de- 
termined a projected separa tion o f l'/2, or 2650 AU, assuming 

~ ; ~ r~„ „ , a distance of 2.2 kpc. Chini et al] J2004I) assumed that the sec- 

Sencl offprint requests to: R. Follert , . , XT ~ . ., . .. . ., , ■ -c , 

* _ , , ond component (No. 2 in their paper) contributes a significant 

Based on observations within the ESO programs 075.C-0755(AB), . £ a , ., , , • £ , , ., 

„__ r „ . .„ . Rr . ^ c amount of flux only in the visual and near infrared wavelength 

** " ., , ., ' T . .. ... n . , „ , c , , c regime. Here, we focus our r esearch on the brighter compo- 

Fellow ot the International Max Planck Research School for ,„, ■, ■ TT^~- — : ri u nA , i „ , „, . . , , 

Astronomy and Cosmic Physics at the University of Heidelberg nent < Na 1 m| Chini et al. || 2004|) . Furthermore, Chini et al. used 

(IMPRS-HD) multi - color photometric data, ranging from the Gunn i band to 
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the MIR, and a TIMMI2 spectrum to fit a simple radiative trans- 
fer (RT) model, assuming a spherical dust distribution. Due to 
excess in the infra r ed pa rt of the spectral energy distribution 
(SED), IChini et al.l d2004 proposed that the KWO is a heav- 
ily embedded Her big Be star, in agreement with the work of 
Port er et al.l (Q~998). In this paper, we present new MIR inter- 
ferometry for the KWO, and we re - evaluate the SED for this 
source, based now on 2D RT simulations. We discuss the statis- 
tics of the model parameters, and how the interferometric results 
constrain these parameters. 



2. Observations and data reduction 

The interferometric data were obtained with the MIDI instru- 
ment installed at the VLTI dLeinert et alj 120031) . The KWO 
was observed during five nights at five different projected 
baseline /projected angle configurations in the N-Band (8- 
13/im). For the first three nights of observation (2005-06-24, 
2005-06-26, 2006-05-18), the 8.2 m diameter Unit Telescopes 
(UTs) were used. For the other two nights, the 1.8 m Auxiliary 
Telescopes (ATs) were used for the observation. At the time of 
submission, this paper represented the first published analysis 
of MYSO data using the ATs. 



the interferometric signal or correlated flux. While measuring 
the correlated flux, the optical path difference (OPD) is first 
changed on a coarser scale (some millimeters) to search for the 
interferometric fringes. Once they are found, the OPD is swept 
continuously around the zero OPD (the white light fringe) with 
an amplitude of approximately 80 //m at a frequency of a few 
hundred scans in three to four minutes. After the correlated flux 
is determined, the photometric signal (or uncorrected flux) is 
measured. 

For that purpose, the light from one telescope is blocked, 
allowing measurement of the flux from the other and vice versa. 
During the recording of the photometric signal, chopping with 
a throw of 10" and a frequency of roughly 2 Hz is applied 
to remove the background signal. Appropriate calibrator stars 
were observed in the same way immediately before and after 
the KWO (see table Q] for details), in order to account for 
instrumental effects which distort the visibility curve. 

The MIDI prism was used as the dispersing element, providing 
a spectral resolution of R w 30. The data obtained with MIDI 
were first reduced using the MIA+EWS package (VI. 5.1), 
developed at the MPIA Heidelberg and the Leiden Observatory 
(see http : //www .mpia-hd .mpg . de/MIDI/ or Jaffe (2004)). 



Table 1. List of observed objects. All objects other than the KWO 
are calibrator stars (see above). The objects marked with an asterisk 
are included in the subsequent analysis. See text for details. The third 
column, labelled B, lists the projected baseline in meters. 



Object 


Date and 
Time 


B 

[m] 


PA[°] 

EofN 


UTs 


HD 107446* 
HD168415* 
HD169916 
HD169916 

KWO* 
HD 1699 16* 
HD177716 


2005-06-24T03:24:09 
2005-06-24T04:22:42 
2005-06-24T07:04:18 
2005-06-24T07:59:14 
2005-06-24T08:23:34 
2005-06-24T08:48:54 
2005-06-24T09:59:36 


61.51 
62.34 
56.89 
52.06 
43.87 
47.75 
47.07 


167.2 
108.6 
127.4 
138.1 
140.7 
150.3 
160.2 


3-4 

3-4 
3-4 
3-4 
3-4 
3-4 
3-4 


HD 1699 16* 
HD 107446 
HD169916 
HD145544 
KWO* 


2005-06-26T00:49:52 
2005-06-26T01:44:57 
2005-06-26T02:29:52 
2005-06-26T03:58:14 
2005-06-26T04:25:ll 


55.84 
38.57 
55.88 
43.41 
55.79 


-6.5 

57.5 
8.0 
44.6 
24.3 


1-2 

1-2 
1-2 
1-2 
1-2 


HD95272* 
HD169916 
HD101666 
HD120404 
HD 120404 
HD168415 

KWO* 
HD 1699 16* 


2006-05-18T00:36:17 
2006-05-18T04:07:17 
2006-05-18T05:41:02 
2006-05- 18T06:42:24 
2006-05-18T07:36:18 
2006-05-18T09:27:57 
2006-05- 18T09:59:02 
2006-05-18T10:28:38 


46.63 
43.99 
26.03 
28.96 
25.98 
44.99 
43.52 
40.01 


43.4 
10.4 
58.0 
100.1 

114.7 
47.8 
48.4 
52.2 


2-3 
2-3 
2-3 
2-3 
2-3 
2-3 
2-3 
2-3 










ATs 


HD150798 

KWO* 
HD211416 
HD168454* 


20()6-08-09T23:30:24 
2006-08-10T01:27:29 
2006-08-10T01:52:18 
2006-08-10T02:12:13 


31.6 
31.7 
31.5 
31.8 


65.6 
69.6 

22.6 
73.9 


DO - GO 
DO -GO 
DO -GO 
DO - GO 


KWO* 
HD211416 
HD150798 
HD150798 
HD150798 
HD206778* 


2006-08-13T00:22:43 
2006-08-13T00:48:05 
2006-08-13T01:40:00 
2006-08-13T02:46:00 
2006-08-13T03:42:34 
2006-08-13T04:35:19 


15.1 
15.7 
14.9 
14.2 
13.6 
15.4 


64.76 
9.9 
-81.7 
-65.5 
-50.7 
74.8 


E0-G0 
E0-G0 
E0-G0 
E0-G0 
E0-G0 
E0-G0 



Table Q] lists the observations, dates, times, telescopes used, as 
well as the associated baseline lengths and orientations. The 
observations followed the procedure described in Leiner tet al.l 
(2004). MIDI was set to work in HighSens - mode. The beams 
from the two telescopes were first combined to determine 



20 ro- 
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Fig. 1. Photometry of the KWO. The dashed, gray lines depict the spec- 
tra as obtained by MIDI, the dash - dotted, dark gray line a rchival IRS 
data, the diamonds the N-band SED from lChini et alj d2004h . and the 
bold, black line the spectrum used for further analysis. The vertical gray 
bar denotes the wavelength regime of the atmospheric ozone absorption 
feature. 



We found that the photometry results, when calibrating with 
different calibrators belonging to one night of observation, were 
not consistent (see Fig. [TJ. This is a clear sign that the nights 
of observation were far from photometric. Hence, we employed 
a different data reduction scheme. First, the correlated fluxes 
of the KWO for each night of observation were determined 
and afterward calibrated by the corresponding set of correlated 
fluxes of the calibrator stars (see Fig. |2|. The correlated flux is 
less affected by atmospheric emission, due to the differential 
recording technique. For determining the correlated flux, two 
subsequent 180° phase - shifted measurements are subtracted 
from each other. Due to the high fringe scan speed, the back- 
ground is subtracted far more efficiently than by the chopping 
procedure. Accordingly, the variations in the correlated fluxes 
were negligible within the errors. The only exception are 
artifacts caused by the atmospheric absorption features (e.g. 
the ozone feature at 9.4-9.9//m). If the science target and the 
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calibrator are observed at different air masses, the calibrator 
will be more (or less) affected by the atmospheric absorption, 
which may result in emission - (absorption -) like artifacts in the 
calibrated correlated fluxes. Hence, only calibrators observed 
at air masses comparable to those of the KWO were taken into 
account for the subsequent analysis (marked with an asterisk in 
table [TJ. Nevertheless, the correlated fluxes in the wavelength 
range covering the ozone feature might still be slightly affected; 
hence, we omitted these data from the analysis. 
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Fig. 2. The measured correlated fluxes. 



12 3 4 5 

Pixel steps in arcseconds 

Fig. 3. Comparison between radial cuts of image profiles of the KWO 
and the calibrator star HD 1699 16, derived from MIDI UT acquisition 
images. 



In addition to the interferometric data, we also took acquisition 
images from the UT measurements at 11.8/im. After back- 
ground subtraction, we compared the profiles of the science 
target and a calibrator star. They are equal within the errors, 
which is remarkable since the observation of the calibrator 
has a clearly better adaptive optics correction (on - source 
guiding). Hence the KWO is not resolved in the MIR with 
single 8m- class telescopes (see Fig. [3}. 



In order to deduce the calibrated visibilities, the calibrated cor- 
related fluxes need to be divided by the calibrated uncorrelated 
flux (i.e. the spectrum). 

Since we wanted to avoid the non - photometric data obtained 
by MIDI, we instead used TIMMI2 data (see Fig. Q] solid 
bold black line). The ESO archive contains mid-infrared grism 
spectroscopy o f the KWO between 8.0 and 13.0 //m obtained 
with TIMMI2 dReimann et alJ l2000h with a nominal spectral 
resolution of 230. The observations were performed within the 
ESO programme 71.C-0185(A) on July 25, 2003. Using a slit 
width of l'/2, the measurements were done with chopping and 
nodding throws of 10", with a total integration time of seconds. 
The star HD 1699 16 served as spectropho tome trie standard 
star. For the re - reduction of the s e data, we used the TIMMI2 
pipeline by ISiebenmor gen et all d2004l) . The same spectrum 
was taken into account for the subsequent SED fitting. For that 
purpose, the data were re - sampled to matc h the wavelength gri d 
used by the Online SED model fitter dRobitaille et all 0007). 
see also section [3~TTi. 



Table 2. Flux densities of the KWO 



Furthermore, S pitzer IRS data fro m the archive and the SED 
as specified in IChini et al.l d2004l) complemented the spectral 
information. The post -basic calibrated IRS data (Request 
Key 11546624, PI M. Wolfire) was reduced with the standard 
software version 15.3.0 and subsequently re - binned. The KWO 
is quite luminous in the NIR. Thus, the Spitzer IRAC frames 
are overexposed. There are no Spitzer MIPS data of the KWO. 
Unfortunately, no (sub-) millimeter continuum data are available 
for the KWO. The SED used in the subsequent analysis is 
summarized in table [2] For the error analysis, variations in the 
correlated flux (during one measurement cycle) as specified 
by the data reduction tool were taken into account, as well as 
the error specified for the TIMMI2 spectrum. If more than one 
calibrator provided usable data, the respective visibility curves 
were finally averaged. 



wavelength 


F v 


°"F„ 


aperture 


data 


[/an] 


[mJy] 


[mJy] 


arcsec 


source 


0.8 


0.025 


0.005 


2 


Chini et al. (2004) 


1.25 


7 


1 


2 


Chini et al. (2004) 


1.65 


91 


9 


2 


Chini et al. (2004) 


2.2 


550 


60 


2 


Chini et al. (2004) 


3.9 


4400 


1000 


2.5 


Chini et al. (2004) 


4.4 


9700 


4500 


18 


MSX B (archive) 


4.6 


5800 


1500 


3 


Chini et al. (2004) 


8.06 


22920 


930 


1.2 


TIMMI2 


8.50 


28890 


1180 


1.2 


TIMMI2 


8.97 


26670 


1090 


1.2 


TIMMI2 


9.44 


24110 


1030 


1.2 


TIMMI2 


9.95 


22310 


940 


1.2 


TIMMI2 


10.49 


25920 


1080 


1.2 


TIMMI2 


11.05 


29510 


1240 


1.2 


TIMMI2 


11.65 


31450 


1320 


1.2 


TIMMI2 


12.27 


31730 


1380 


1.2 


TIMMI2 


12.93 


31090 


1370 


1.2 


TIMMI2 


14.7 


33070 


1650 


18 


Chini et al. (2004) 


17.8 


55488 


13900 


5 


Chini et al. (2004) 


21.3 


68760 


3440 


18 


Chini et al. (2004) 


23.04 


105000 


6100 


12 


IRS 


26.96 


129000 


5900 


12 


IRS 


31.56 


156000 


7100 


12 


IRS 



3. Results 

The analysis of spectrally dispersed visibility curves is non- 
trivial. Analytical solutions can only be found if one assumes 
very simple models and/or geometrical shapes. For this reason, 
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we first fit the SED to pre - calculated models. The set of param- 
eters which describe the respective model is then used as input 
for a Monte Carlo radiative transfer model. The resulting syn- 
thetic images are finally Fourier - transformed in order to pro- 
duce synthetic visibility curves which can be compared to the 
measurements (Lin z et al.ll2009t) . 

3.1. SED fitting 

The SED data served as input for the Online SED model fitte£] 
(Robitaille et aD l2007l) . The fitter is based on a grid of 200,000 
pre - calculated, 2D Monte C arlo radiative t ransfer models com- 
puted with the code of Whitn ey et al.l (l2003). Each of these mod- 
els is characterized by a set of 28 param eters, which are not all 
independent of each other (for details see Robitaille et al. (2006) 
and Robitaille (2008)). The models include the properties of the 
central object (mass, radius, temperature), the envelope (accre- 
tion rate, outer radius, inner cavity density and opening angle), 
the disk (mass, inner and outer radius, accretion rate, scale height 
factor, flaring angle) and the ambient density. These parameters 
have been varied to cover a large range of possible configura- 
tions of young stellar objects, yet the parameter volume is not 
covered uniformly. The result is a first grid - like set of 20,000 
models. For each of these models, the SED was calculated for 
10 different inclination angles (defined between the axis of sym- 
metry of the model and the line of sight), distributed equally over 
the cosines of the inclination angle. Hence, the final grid consists 
of 200,000 models. The grain model used consists of a mixture 
of astronomical silicates and graphite i n solar abundance, w here 
the optical constants are taken from lLaor & Draine ( 1993J). In 
addition, the optica l properties are averaged over the size distri- 
bution (taken from Kim et al. 1994) and composition. For further 
details see Robitaill e~t alj (120061) . 

In addition to the SED, the user is required to specify some as- 
sumptions on the distance and interstellar extinction toward the 
object as input for the Online SED model fitter. These assump- 
tions are treated as boundary values in the fitting procedure. The 
fitter then performs a global x 1 minimization procedure to eval- 
uate which parameter combination fits the observed SED best. 
The results are presented ordered by increasing deviation be- 
tween the modeled and specified SEDs. For each of these mod- 
els, a parameter file is available which contains the set of param- 
eters describing this particular model. These files then serve as 
input to the Monte Carlo radiative transfer code. 

For this work, the distance and the inte r stellar extinction were 
chosen in accordance with IChini et all d2004l) . who proposed 
a distance d = 2.2 kpc and a total visible extinction of A v at 
24 mag, including extinction due to both the intervening in- 
terstellar matter and the l ocal circumstellar struc tures. The re- 
cent re - determination by iHoffmeister et al.l d2008l) . who found 
d = 2.1 +0.2 kpc for the M17 SW star cluster, was also taken 
into account. Hence, we adopted as boundaries for the fitting 
procedure d = 1.9 - 2.3 kpc, while the interstellar extinction was 
assumed to be larger than 5 mag. We want to point out again 
that the distance, the interstellar as well as the total extinction 
are treated as free parameters (within the boundaries). Figure [4] 
presents the ten best SED fits. 

The Spitzer IRS data were considered as upper limits. There will 
be a significant contribution to the luminosity by the interstellar 
material surrounding the KWO, due to the much larger aperture 
of IRS at these wavelengths. Apparently, the fitted curve overes- 



timates the flux at 3.9 and 4.6 fim. In order to assess the result, we 
ran the SED fitter again, including each data point shortward of 
8 /im twice. The resulting set of solutions indeed involves some 
spectra which agree on the relatively low flux at 3.9 and 4.6//m. 
Yet, we found that all these solutions cannot account for the pro- 
nounced silicate feature (see Fig. [T). In addition, we included a 
MSX B2 data point from the archive. Due to the uncertain cali- 
bration of the MSX B2 data, the associated error is large; hence, 
it does not have a strong effect on the SED fit. Nevertheless the 
MSX B2 flux and the best fit match quite well. Hence, we kept 
the result as presented in Fig. [4] 



<o 10" 



.10" 
10" 
10" 
10" 





_ 


iii! 




' 






1 1 


■ 




, , , 1 


■ 
















KWO 








1 


-7 






















1 


-8 








/•» 






•m%* # 


WTW 






1 


-9 


















\s 




1 


10 






















-. 


11 






1/ 


\ 














\ I 












\ 














■■ 


12 




/ 
/ 








\ 


\ 










f 


13 
























1" 


14 














. i 




\ 


, , , 1 







10 
A (/im) 



100 



http : //caravan . astro . wise . edu/protostars/index . php 



Fig. 4. Best SED fits, obtained with the Online SED model fitter. The 
black solid line corresponds to the best fitting model, the gray lines to 
the 9 next best fits. The dashed line shows the stellar photosphere model 
with the best - fit interstellar extinction taken into account. The dots and 
triangles represent the data the fit was applied to, dots denote regular 
data points and triangles upper limits, respectively. 



The compiled SED used for this work is based on data cover- 
ing the peak and Wien side of the SED. No published data ex- 
ist for the Rayleigh-Jeans part, which could distinguish KWO 
from the surrounding high level of extended (sub-) millimeter 
flux present in the M17 SW region. Therefore, the SED of KWO 
is not constrained at far infrared and submm wavelengths, and 
several model configurations fit the observed optical /infrared 
SED almost equally well. For this reason, a statistical approach 
to constrain the meaningful parameter space was adopted (see 
Sec. @). We created histograms of the different model param- 
eters (accounting for the 1000 best fitting models), with each 
count weighted by l/x 2 - High densities in the respective param- 
eter space hence give an indication of the probability that a cer- 
tain parameter will give a value in the associated range. 

3.2. Radiative transfer modeling and synthetic visibilities 

For the 10 best -fitting configurations, MCRT simulations were 
performed using the ttsre code of Whitnev et al. (2003). The 
simulation was run each time using a total of 100 million test 
photons. We modified the MCRT code for our purposes by defin- 
ing 1 1 narrow filters centered at equally spaced wavelengths cov- 
ering the N - band. Hence, we obtained for each model 1 1 syn- 
thetic images. Figure [5] shows an example of the frames created 
by this procedure. We combined 3 synthetic images of different 
wavelength to create an RGB composite. In addition, we applied 



R. Follert et al.: Mid - infrared interferometry of massive young stellar objects II 



a logarithmic intensity scaling in order to reveal the diffuse outer 
regions of the structure. The picture was adaptively smoothed to 
reduce the graininess inherent to Monte Carlo simulations. 

Each of these synthetic images was first embedded in a larger, 
empty frame, since the MCRT code only allows for a certain spa- 
tial resolution in combination with a respective maximal spatial 
outer scale. Then, each image was Fourier- transformed in order 
to obtain the related (u , v) - spatial frequency spectrum. Due to 
the increased spatial outer scale, the resolution in Fourier- space 
was refined, which in turn helped to avoid numerical discretiza- 
tion effects. 



clearly resolved on all baselines. Moreover, the silicate absorp- 
tion feature we see in the total intensity spectra is not seen in any 
of the visibility curves. This is due to the fact that the feature is 
equally strong in both the correlated and uncorrected flux. Since 
we see the same effect in every visibility curve, although they 
are measured on clearly distinct baseline lengths, this indicates 
that the absorbing dust is situated in foreground structures not 
directly associated with the KWO. The existence of such a dust 
screen in f ront of the KWO w as already proposed by Chini et all 
(l2004t) and lPorter et alJd 19981) . who suggest that the KWO might 
actually be located behind the SW obscuration of M17. 




Fig. 5. Synthetic RGB composite image for model 3003929 (see Table 
[5J at an inclination angle of 75°. The 12.5 /im bin is depicted in red, the 
10.5 //m bin in green and the 8.5 /im bin in blue. The total visible area 
isl200x!200AU. 



Next, the Cartesian coordinates of the (u , v) - plane were trans- 
formed to polar coordinates. The synthetic visibility curves were 
constructed by adopting the values from the synthetic (u , v) - 
spectra from each of the respective wavelength bin images at the 
spatial frequencies corresponding to the projected baselines used 
in the observations. There is one remaining degree of freedom: 
while the inclination angle i is fixed due to the SED fitting and 
the respective MCRT modeling, the symmetry axis of the syn- 
thetic images runs north-south, which is not necessarily true for 
the science target. Hence, an additional angular parameter is re- 
quired to account for the observed position angle. This parameter 
was derived by minimizing the deviation between the synthetic 
visibility curves and the measured visibility curves by rotating 
the u - v plane. Hence, we obtained for each model spectrally 
dispersed, synthetic visibility curves at the projected baselines 
used for the observation. These will be compared to the mea- 
sured visibilities in 14.41 



This approach is similar to lde Wit et al.l ( 120 10), who treated the 
MYSO W33A. They use the RT models of lWhitnev etaD (Eool 
to gen erate synthetic visibility curves as well. Yet. Ide Wit et al.l 

(2010) were able to constrain various model parameters in ad- 
vance due to better availability of data for W33A, especially in 
the submm wavelength range. 

3.3. Interferometric observations 

The measured calibrated dispersed visibilities appear in Fig. 
[6] The depicted errors represent 3<x. The visibility amplitudes 
are smaller than unity for all observations. Hence, the KWO is 




x W 



Fig. 6. Calibrated visibility curves. 



4. Discussion 

4.1. Statistical reduction of the parameter volume 

SED data fitting of MYSOs is in general ambiguous. In order 
to investigate how the observational errors translate into varia- 
tions of th e RT model parameters, we followed the approach of 
Pinte et al. (2008). Note that the histograms (see section [3~TT ) rep- 
resent the distribution of the respective parameters taken from 
the set of the 1000 best -fitting models as obtained from the 
Online SED model fitter, no a priori assumptions were involved. 



4.1 .1 . The stellar properties 

Figures [7] and [8] show the distribution of the fitted stellar masses 
and temperatures, respectively. The fitted stellar temperature is 
strongly peaked around 32000 K. The stellar mass distribution is 
peaked around 15 M , with a side lobe towards 10 M Q . Hence, 
the most probable embedded star for the KWO would be an early 
B star. 



4.1.2. Disk and envelope 

Figure|9]plots the envelope accretion rate. It is interesting to note 
that most of the models fitting the SED have no envelope at 
all: these models consist of the central star and a disk only. In 
the Robitaille et al. (2007) models, this is denoted by setting the 
envelope accretion rate formally to zero. This is in contrast to 
IChini et all (|2004 who fit the SED with a spherical dust distri- 
bution with a radial density profile varying as r , which might 
be interpreted as an envelope - like structure. Due to the inabil- 
ity of this approach to properly reproduce the silicate feature, 
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Fig. 7. Distribution of the stellar mass [M Q ] 



Fig. 10. Distribution of the disk mass [M e 




2.0x10* 2.5x10* 3.0x10* 3.5x10* 4.0x10* 

Stellar Temperature [K] 

Fig. 8. Distribution of the stellar temperature [K] 

one of their suggestions was a disk - like structure, which would 
produce silicate features both in emi ssion a nd absorption, de- 
pending on the disk orientation. Chini et all (120041) identify the 
reflection nebula surrounding the KWO with this spherical dust 
distribution. The arguments made above and our interferometric 
measurements (see sections 14.31 and I4.4t favor a configuration 
without a small - scale optically thick dust shell directly envelop- 
ing the KWO. Thus, the scattering material in the reflection neb- 
ula is probably not bound to the KWO itself. Instead, we suggest 
that the reflecting material is either a remnant of the natal molec- 
ular cloud or it is situated in front of the KWO a long the line of 
sight (as already proposed by Port er et al.l (Il998l) . 



no 10"" 10"'° 10" s 10"" 10" 7 10" B 10~ 5 10"" 1C 
accretion Envelope Accretion rate [M e /yr] 

Fig. 9. Distribution of the envelope accretion rate [— ] 



Figure \TU\ shows the histogram of disk masses. We find 



that configurations with a disk of intermediate mass (10~ 3 to 
1CT 1 M ) are most probable. Nevertheless, we want to stress that 
(sub-) mm data with sufficient angular resolution are not avail- 
able. Only such observations would give the total mass of cir- 
cumstellar material by measuring the optically thin thermal dust 
emission. Without such data, any assumptions on the mass of the 
circumstellar material will be model dependent. 




9 11 13 15 '7 19 21 23 25 
nner Disk Radius [R,„ B ] 



Fig. 11. Distribution of the Disk inner radius [R su b] model parameter. 
Values of - 1 correspond to envelope - only models 



Figure \TT\ shows the corresponding inner radii of the proposed 
disk structure, in units of the thermal dust sublimation radius. It 
is interesting to note that the most probable inner radii are larger 
than the dust destruction radius. This result is roughly consis- 
tent with the relation b etween the total lum i nosity and the source 
size as established by Millan-Gabet et al. (2007). However, the 
radius of the best - fitting model (Sect. 14.41 table [3]) exceeds the 
dust sublimation radius by a factor of two, suggesting ongoing 
inner disk dispersal by photo - evaporation and/or stellar winds. 



4.2. Object models from the SED Fitting 

As a result of the SED fitting, we obtained 4 slightly different 
model configurations with various inclination angles which re- 
semble the measured SED best (see table|2]). In table[3j we sum- 
marized the most important structural parameters. A more de- 
tailed overview can be found in table[5]in the appendix. 

All of these model configurations consider a disk only; hence, 
no additional circumstellar material is needed to account for the 
near and mid - infrared fluxes. The central object is in every case 
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quite massive. According to iMartins et al.l (120051) . the objects 
can be classified as early B stars. The disks themselves are mas- 
sive as well, yet the distance of the inner rim to the central object 
differs considerably. Since the inner radius is directly related to 
the spatial distribution of the dust and, therefore, to the spatial 
intensity distribution, we should be able to distinguish the dif- 
ferent models when comparing the synthetic and measured visi- 
bility curves. We explicitly do not specify an outer scale for the 
disk, since our data are not capable of confining the size of the 
disk in a model - independent way. 



Table 3. Parameters and derived quantities of the models. Note that 
these are disk -only models. 







Grid Model Number 






3003929 


3005903 


3008813 


3012934 


stellar mass [M Q ] 


15 


16 


16 


14 


stellar temperature 
[K] 


31300 


32300 


32200 


30500 


disk mass [M Q ] 


9xl0~ 2 


4 xlO" 2 


2xl0- 2 


6 xlO- 2 


disk inner radius 


26 (2.2) 


84 (6.4) 


99 (7.7) 


55(5.1) 


disk accretion rate 

[M /yr] 

Luminosity 

[L Q ] 


7.8 xlO- 8 


3.0 xlO" 6 


2.2 xlO" 7 


4.9 xlO" 7 


20.7 xlO 3 


25.9 xlO 3 


25.0 xlO 3 


17.7 xlO 3 



4.3. General discussion of the visibilities 

Before actually comparing the synthetic to the measured visibil- 
ity curves, we will discuss a few more general conclusions which 
can be drawn from the overall shape of the measured curves 
themselves. The most intriguing fact is that although the base- 
line lengths for the 2005-06-24 and 2006-05-18 measurements 
are almost equal, the respective visibility amplitudes are clearly 
distinct, even within the depicted 3 cr errors (see Fig. [6}. The re- 
spective projected angles (141° and 48°) differ by approximately 
90°. This is a clear indication for a non - rotationally symmet- 
ric intensity distribution, i.e. the mid -infrared emission is dom- 
inated by a flattened structure, e.g. a disk. This interpretation of 
the visibility curves is consistent with the disk -only scenario 
proposed by the statistical analysis of the SED fitting. 

To give an estimate of the size of the source, we fitted a Gaussian 
intensity profilq^ to the visibility amplitudes at 8, 10.5 and 
13yum (see Fig. Q~2] fit at /l=10.5/im). The corresponding full 
width at half maximum were determined as 33, 38 and 40 mas 
(corresponding to 69, 80 and 84 AU, if a distance of 2.1 kpc 
is assumed), respectively. Note, though, that since obviously a 
Gaussian profile is a bad fit to the data, these numbers should be 
rather regarded as rough estimates. 

For instance, if one assumes an extended source of 8 - 1 3 yum ra- 
diation, such as a halo or hot dust in an ultra compact Hn re- 
gion, one would expect to measure rather small visibility ampli- 
tudes. These sources typically have sizes larger than 5xl0 15 cm 
(half -intensity radius, ( Church well et al. 1990), corresponding 
to ss 160 mas at 2.1 kpc), which would imply much lower visi- 
bilities at these baselines due to over- resolving the target (the re- 
spective FWHM in Fig.[T2]would be narrower by approximately 
a factor of 4 to 5). This is in accordance with the results found 
in the previous section, namely that the most probable solution 
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Fig. 12. Fit of a ID Gaussian intensity distribution to the visibility am- 
plitudes at 10.5 yum. 



for the SED fitting is an evolved, disk - only system without a 
dust envelope. The argumentation holds as well for every source 
of N - band emission being much larger than the scales found 
above. It applies, for example, also to a system consisting only of 
a central star plus flared circumstellar disk, seen edge - on. There, 
the direct view onto the hot inner rim of the disk is blocked, and 
we just receive more diluted emission on larger scales. Given our 
assumption that the projected intensity distribution of the KWO 
is not rotationally symmetric, which excludes a face - on disk, 
we would expect some intermediate inclination angle. Indeed, 
the ten best-fitting models are inclined by 57 to 81° (where the 
resulting configuration still allows a direct line of sight to parts 
of the central disk edges (see Fig. [5}. 

4.4. Comparison between modeled and measured visibility 
curves 

Here, we discuss the comparison between the modeled syn- 
thetic visibility curves and the measured ones (see Figs. [L3l and 
fT4l . Apparently, there is quite a discrepancy between the syn- 
thetic and the measured spectrally dispersed visibility curves. 
Nevertheless, there are clear differences in the qualitative shape 
of the modeled visibility curves, which enable us to exclude 
some of the models which otherwise fit the SED of the KWO 
almost equally well. Hence, this approach offers the ability to 
reduce the ambiguities arising from SED - only fitting. 

The difference between measurement and model can have sev- 
eral origins. First, the flux measurements used with the Online 
SED model fitter were incomplete. Lacking submm/mm data, 
we were not able to constrain the flux from the cold dust in the 
outer parts of the disk. This relates directly to the ambiguity in 
the results of the Online SED model fitter (see Fig. [4) The fitted 
SEDs differ considerably in this wavelength region). In addition, 
the modeling still suffers from many imperfections. For example, 
the model u s ed her e does not account f or a p uffed - up inner rim 
dNatta et all (l200lh . iDullemond et all d200lln . yet most of the 
flux in the N - band stems from regions of the disk which would 
be affected in dimension and shape if this were implemented. A 
change of the geometric shape of this region will affect the shape 
of the synthetic visibility curves. The dust model used by the RT 
simulation code is also not well adapted to this disk - only sit- 
uation, since it assumes a dust composition that resembles the 
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typical composition of the halo, not including bigger dust parti- 
cles which are likely present in an evolved disk. Additionally, a 
bigger variety of dust particles will lead to a smoothing of the 
dust destruction zone. Smaller dust grains in this region reach 
their evaporation temperature at larger distances from the cen- 
tral object than bigger ones. Hence, the inner radius of the disk 
as a function of height correlates with the size of the dust grain 
constituents, leading to a more curved, rather than fla t, hot inner 
rim dTannirkulam et all (120071) . Ilsella & Nattal (120051) ). 

Proper error handling is difficult to achieve when comparing 
the synthetic curves with the measurements. A conventional ap- 
proach would be a global x 1 measure, including deviations on 
all 5 baselines. We found that many models fit all the measure- 
ments except the 55.8 m baseline, where they fail to reproduce 
the rather low visibility amplitude. A conventional x 2 measure, 
on the other hand, would overweight these measurements due to 
their small absolute errors. Yet, we found that the relative error 
of all measurements is almost constant. Hence, instead of the ab- 
solute error, we used the square of the relative error as a weight 
in the error analysis. As a result, each curve is almost equally 
significant. In order to justify our ansatz, we did the same analy- 
sis excluding the 55.8 m baseline. We found that in this case, the 
two ways of handling the error deliver the same results, while 
they deviate partially when we included all measurements in 
the standard deviation estimation. In this case, the absolute error 
weighted result was biased towards the 55.8 m baseline. Hence, 
we decided to use the square of the relative error as the weight 
for comparing our measurements to the model. The results of 
this analysis are summarized in table|4] For a better comparison, 
the inclination angle and the x 1 obtained from the SED fit are 
specified for each model in the first two columns. In column 4, 
the position angles (see section l3~2l > best -fitting the data and in 
column 5 the respective estimates of the deviation between the 
synthetic models and the data are given. 



Table 4. Standard deviations of the SED fits and estimates of the devi- 
ation between modeled and measured visibility curves 



Grid model 


Inclin. 


SED Fit 


Pos. 


Visibility 


number 


angle [°] 


x 2 


angle [°] 


9 

X" 


3003929 


76 


141 


16 


73 


3005903 


81 


189 


167 


248 




63 


190 


176 


125 


3008813 


70 
76 


153 

123 


173 
171 


145 
155 




81 


151 


171 


210 




57 


194 


84 


260 


3012934 


63 
70 


199 
189 


89 
92 


284 
317 




76 


188 


92 


338 



The best agreement between model and data is obtained with 
model 3003929, the deviation estimate being almost a factor 
of two smaller than the next best one. It is interesting to note, 
that model 3003929 is not the best SED fit, although the respec- 
tive values of the standard deviation of the SED fit for model 
3003929 and the best fit are comparable. Model 3005903 and 
the four different variations of model 3012934 agree least with 
the data. Hence, we exclude them from the further discussion. 
Note that the four realizations of model 30088 1 3 differ mainly in 
the inclination angle. The comparison between these four mod- 
els and the data shows a trend toward smaller inclination angles. 
Yet, due to the small differences between the four realizations, 



the ones with inclination angles of 63°, 70°and 76° can be con- 
sidered to agree equally well with the data. A comprehensive 
overview of the parameters determining model 3003929 (as well 
as the other models) can be found in Table |5]in the appendix. 

Although model 3003929 agrees best with the data, there are 
obvious deviations between the synthetic and the measured visi- 
bility curves. To begin with, the synthetic visibility amplitude for 
the 55.8 m baseline is higher than the measured one. A possible 
reason for this was given above: we consider this deviation to 
be systematic. There is good agreement for the other baselines, 
yet we find a qualitative discrepancy. Model 3003929 shows in- 
creased visibility amplitudes for the 8 to 9 /mi regime. We think 
that the source of this behavior is directly related to how the 
RT code handles the inner rim. This rather narrow, hot region, 
from which emission in this wavelength regime emanates, is not 
as resolved in comparison to the other models due to the small 
inner disk radius of 3003929. This leads to higher visibility am- 
plitudes. If one assumes that, in reality, a puffed -up inner rim 
exists, the spatial region of the disk emitting in the spectral range 
between 8 and 9 /mi might be considerably more extended, lead- 
ing to the lower visibility amplitudes we see in our measure- 
ments. 



5. Conclusions 

In this paper, we present VLTI observations of the Kleinmann - 
Wright object with baselines ranging between 15 m and 56 m. 
Our target has been resolved for all baselines, with average vis- 
ibility amplitudes between 0.01 and 0.69. In order to analyze 
the interferom etric data, we perfo rmed SED fitting using the 
web tool of Robitaille e t al.l (120071) . The resulting information, 
together with a Monte Carlo radia tive transfer simulation code 
developed bv lWhitnev et al.l (12003). was then used to create syn- 
thetic visibility curves. We compared these synthetic curves to 
the measured visibilities, allowing us to substantially narrow 
down the range of suitable solutions obtained by mere SED fit- 
ting. 

While massive YSOs are usually thought to be associated with a 
significant circumstellar envelope which often dominates their 
infrared appearance, KWO does not show indications for a 
strong envelope component, and hence is probably more ad- 
vanced in its ev olution t owards the main sequence. This is in 
accordance with Porter et al. (1998), who proposed that the clas- 
sification as an embedded MYSO might be misleading, since it 
could be behind M17 SW but not actually embedded in it. We 
found that the silicate absorption feature is missing in all visi- 
bility curves. This indicates that the silicate absorption feature, 
which is clearly visible in the total flux spectrum, is mainly due 
to foreground structures not associated with the KWO. In this 
geometry the correlated and uncorrelated fluxes are equally af- 
fected; hence, the visibilities don't show any pronounced silicate 
features. In addition, our statistical interpretation of the SED fit- 
ting also indicates the absence of a dusty envelope. We find that 
the visibility curves in general can only be interpreted using an 
intensity distribution deviating from spherical symmetry, which 
hints at a circumstellar - disk scenario. Finally, we deduce that a 
flared disk, inclined at 76° and with a mass of »0. 1 M and an 
inner rim radius of 26 AU, reproduces the combination of SED 
fitting and interferometric data best. To obtain better constraints 
for the dimensions of the disk and the amount of circumstellar 
material in general, high - resolution submm/mm data will be 
needed. 
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Fig. 13. 

3cr. 



Comparison of the observed and synthetic visibilities; each figure represents one baseline configuration. The error bars are 
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Table 5. Summary of model parameters. For details on the pa- 
rameter definitions see Robitaille et al. (20Qa). 





3003929 


3005903 


3008813 


3012934 


Inclination 
Angle[°] 


76 


81 


63 


70 


76 


81 


57 


63 


70 


76 


SED fitting 

2 

x~ 


141 


189 


190 


153 


123 


151 


194 


199 


189 


188 


distance 
[pc] 


1905 


1995 


2291 


2291 


2291 


2188 


2188 


2188 


2089 


1995 


Ay 

(total) 


134 


587 


19 


18 


17 


94 


19 


18 


19 


27 


Ay 

(interstellar) 


13 


12 


19 


18 


17 


12 


19 


18 


19 


16 


stellar 
radius [R Q ] 


4.9 


5.1 


5.1 


5.1 


5.1 


5.1 


4.8 


4.8 


4.8 


4.8 


stellar 
mass [M ] 


15 


16 


16 


16 


16 


16 


14 


14 


14 


14 


stellar 
temperature [K] 


31000 


32000 


32000 


32000 


32000 


32000 


31000 


31000 


31000 


31000 


Luminosity 
(total, rLnl) 


20700 


25900 


25000 


25000 


25000 


25000 


17700 


17700 


17700 


17700 


envelope accr. 
rate [M /yr] 
































disk mass 
[M Q ] 


9x1 0~ 2 


4x1 0~ 2 


2xl0- 2 


2xl0- 2 


2xl0~ 2 


2x1 0~ 2 


6xl0 2 


6xl0 2 


6x1 0~ 2 


6x1 0~ 2 


disk inner 
radius [R„,i] 


2 


6 


8 


8 


8 


8 


5 


5 


5 


5 


disk inner 
radius [AU] 


26 


84 


99 


99 


99 


99 


55 


55 


55 


55 


disk accr. 
rate [M /yr] 


8x1 0-° 8 


3xl0- 6 


2x1 0~ 7 


2xl0~ 7 


2xl0~ 7 


2x1 0~ 7 


5xl0- 7 


5xl0- 7 


5xl0- 7 


5xl0- 7 


disk flaring 
parameter /? 


1.15 


1.04 


1.10 


1.10 


1.10 


1.10 


1.10 


1.10 


1.10 


1.10 



Note that the total extinction given in the table is the sum of the 
interstellar and intrinsic extinction. The intrinsic extinction for 
disk - like structures is given for the direct line of sight; hence, it 
depends considerably on the inclination and flaring of the disk. 
If the total and interstellar extinction are equal on the other hand, 
this indicates a direct line of sight to the central stellar source. 
Robitaill e et al.1 (|2006) state that the disk structure is given by: 



p(m, z) = po 






lrzi2\ 



(zM-m 



where h re xrfi and a — j3 + 1 . 



